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Domain- wall fermions with U(l) dynamical gauge fields 
in (4+l)-dimensions * 

Sinya Aoki and Kei-ichi Nagai a 

a Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 

We carry out a numerical simulation of a domain- wall model in (4+1) dimensions, in the presence of a 
quenched U(l) dynamical gauge field only in an extra dimension, corresponding to the weak coupling limit 
of a (4-dimensional) physical gauge coupling. Our numerical data suggest that the zero mode seems absent in the 
symmetric phase, so that it is difficult to construct a lattice chiral gauge theory in the continuum limit. 



1. Introduction 

Construction of chiral gauge theory is one of 
the long standing problems of lattice field the- 
ories. Because of the fermion doubling prob- 
lem, the lattice field theory discretized naively 
becomes non-chiral. Several approaches on lat- 
tice have been proposed to overcome this difficul- 
ties, but none of them have been proven to be 
successful. 

Recently Kaplan has proposed a domain-wall 
model in order to construct lattice chiral gauge 
theories Q]. The model consists of Wilson fermion 
action in (2n+l) dimensions with a fermion mass 
term being the shape of a domain wall in the 
(extra) (2n+l)th dimension. In the case of free 
fermions it is shown for < mo < I, where mo 
denotes the domain wall mass height, that a mass- 
less chiral state arises as a zero mode bound to 
the 2n-dimensional domain wall while all doublers 
have large masses of the lattice cutoff scale. 

In a way of introducing dynamical gauge fields 
two variants of this model have been proposed: 
the waveguide model|| and the overlap for- 
mula||. However, it has been reported that the 
chiral zero mode disappears for these two vari- 
ants even in the weak coupling limit, due to the 
roughness of gauge field Q. In the original model 
the roughness of the gauge field is replaced with 
the dynamical gauge field in the extra dimension, 
and in the weak coupling limit of the 2n dimen- 
sional coupling, 2n dimensional links U fl (x.s) = 



1 (/z = 1, ■ • ■ , 2n) while only extra dimensional 
links C/an+i(x, s) become dynamical. We would 
like to know whether the fermionic zero modes 
exist on the domain wall or not in this limit. This 
question has already been investigated in (2+1) 
dimensions B, but the result, in particular, in the 
symmetric phase is not conclusive, due to pecu- 
liarity of the phase transition in the 2 dimensional 
U(l) spin model. Therefore we numerically inves- 
tigate this model with U(l) gauge field in (4+1) 
dimensions and report the result here. 

2. Mean-field analysis 

Before numerical investigation we estimate the 
effect of the dynamical gauge field in the extra 
dimension using the mean-field analysis. In our 
mean-field analysis all link variables in the extra 
dimension are replaced with (x, s)-independcnt 
constant z (0 < z < 1), so that the fermion prop- 
agator is easily obtained [Q ^ 



G(p)s 



iJ2^+ M ( z ) G R (p)P L 



M / 



where p^ = sin(p M ) and P R/L = (1 ± 7s)/2. 
Corresponding fermion masses are obtained from 
Gr and Gl in p — > limit, and we found that 
fermionic zero modes exist only for 1 — z < mo < 
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1. Thus the critical value of domain wall mass is 
m§ = 1 — z and therefore no zero mode survives 
if z = 0. 

3. Numerical analysis 

We investigated the (4+1) dimensional U(l) 
model numerically, using a quenched approxima- 
tion. At the zero physical gauge coupling the 
gauge field action of the model is reduced to the 
2n dimensional spin model with many copies: 



E 



ReTr 



U D (x,8)W D (x + ii,a) , (1) 



where D — In + 1. Therefore there exists a 
phase transition which separates a broken phase 
from a symmetric phase. We calculated the or- 
der parameter v using a rotational technique Q 
and found j3 s = 0.29 corresponds to the symmet- 
ric phase and f3 s — 0.5 to the broken phase. We 
calculated the fermion propagator over 50 config- 
urations at (3 S — 0.5 and 0.29 on L 3 x 32 x 16 
lattices with L —4, 6, 8. For the fermion field 
we take periodic boundary conditions except the 
anti-periodic boundary condition in the 4th di- 
rection. The errors are estimated by the single 
elimination jack-knife method. 

At s = we obtained the inverse propagator 
G^ 1 and G7, 1 for pi,p2,P3 = as a function of 
P4. We obtain fermion mass squared, rnj, by ex- 
trapolating G^ 1 and G^ 1 linearly to p± = 0. 

3.1. Broken phase 

In Fig^ we plotted m/ in the broken phase 
{(3 S = 0.5) as a function of mo- As seen from 
this figure, the finite size effect is small, and the 
left-handed modes are always massive, while the 
right-handed modes are massless if mo is larger 
than about 0.6. Therefore, we conclude that chi- 
ral zero modes can exist in the broken phase. 

3.2. Symmetric phase 

Let us show the fermion mass in the symmet- 
ric phase (p s — 0.29) in Fig.||. In the smallest 
lattice size the chiral zero modes seem to exist. 
However, for large lattices, the mass difference be- 
tween the left- and right-handed modes becomes 
smaller. This suggest that the fermion spectrum 
becomes vector-like in the infinite volume limit. 
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Figure 1. m/ vs. toq i n the broken phase 
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Figure 2. m/ vs. mo in the symmetric phase 

However, since the fermion mass near mo = 1.0 is 
so small, from this data alone, we cannot exclude 
a possibility that the critical mass mg is very close 
to 1.0. 

To make a definite conclusion on the absence of 
chiral zero modes in the symmetric phase, we try 
to fit the fermion propagator using the form of the 
mean-field propagator with the fitting parameter 
z. We show the quality of the fit in Fig.||. This 
figure shows that the fermion propagator is well 
described by the mean-field propagator. 

In Fig.|| and Fig.|| we plotted the parameter z 
obtained the above fit as a function of 1/L The 
parameter z's are almost independent of mo at 
the each 1 /L except for the right-handed ones at 
mo = 0.99. The solid circles represent the order 
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Fi gure 3. G ft vs. sin (p) in the symmetric phase 
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Figure 4. z (right-handed) vs. 1/L in the sym- 
metric phase 

parameter v. The behaviors of z at different mo 
are almost identical each other and are very sim- 
ilar to that of v except the right-handed ones at 
mo = 0.99. This suggest that z can be identified 
with v and therefore z becomes zero as the lat- 
tice size goes to the infinity. If this is the case the 
fermion spectrum of this model becomes vector- 
like in the symmetric phase. 

4. Conclusions 

We have carried out the numerical simulations 
of the U(l) original domain-wall model in (4+1) 
dimensions in the weak coupling limit of the 4- 
dimensional coupling. 

In the broken phase, there exist chiral zero 
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Figure 5. z (left-handed) vs. 1/L in the symmet- 
ric phase 

modes on the domain wall for mo > TOq. The ex- 
istence of the critical mass to§ is predicted by the 
mean-field analysis. On the other hand, in the 
symmetric phase, the analysis using the mean- 
field propagator suggests that this model becomes 
vector-like. We should note, however, that the 
right-handed modes at too = 0.99 behaves differ- 
ently and the similar behavior was also found in 
the (2+1) dimensional modelQ. Therefore in the 
future, we must investigate this point in detail, 
for example, increasing the statistics. 

Besides this point the results from both phases 
suggest that this model becomes vector-like in the 
continuum limit, which should be taken at the 
critical point (: the point between two phases). 
Therefore, it seems also difficult to construct a 
chiral gauge theory on lattice via the original 
domain- wall model. 
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